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1 ' 

We study statistical signatures of composite bosons made of two fermions using 

a ■ 

^ ' a new many-body approach. Extending number-states to composite bosons, two- 

cr 

particle correlations as well as the dispersion of the probability distribution are 
analyzed. We show that the particle composite nature reduces the anti-bunching 
' effect predicted for elementary bosons. Furthermore, the probability distribution 

£\\ • exhibits a dispersion which is greater for composite bosons than for elementary 



bosons. This dispersion corresponds to the one of sub-Poissonian processes, as for a 



CO 

o 

, quantum state, but, unlike its elementary boson counterpart, it is not minimum. In 



general, our work shows that it is necessary to take into account the Pauli exclusion 
j_j ■ principle which takes place between fermionic components of composite bosons 

a 

along the line here used - to possibly extract statistical properties in a precise way. 
PACS numbers: 71.35.-y, 03.75.-b 

Quantum particles obey statistical laws that are highlighted by high-order correlations. 
This has first been shown for light fields, by the discovery of photon bunching in two-photon 
correlations from an incoherent source [l|. Thereafter, second-order correlations have be- 
come increasingly important, particularly to distinguish classical from quantum states of 

light a a. 

Indeed, classical states are characterized by bunched two-photon correlations 
while quantum states yield anti-bunching. For elementary bosons, two-particle correlations 
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are directly connected to the variance of the probability distribution which is usually ex- 
pressed through the so-called Mandel parameter Q (4]. The latter compares the dispersion 
of the distribution to the one of Poissonian processes. Hence, quantum states correspond to 
negative values of Q, i.e., they follow sub-Poissonian statistics, while classical states lead to 
positive values of Q since these can solely follow Poissonian and super- Poissonian statistics. 
Recent advances with ultra-cold atoms have shown that analysis of second-order correla- 

n 

tions is also well suited to probe matter states (see [5| and references therein). Most notably, 
anti-bunching and bunching effects have been observed while measuring two-particle corre- 
lations of degenerate Fermi and Bose gases |f3, B, S|- In quantum mechanical terms, these 
effects are due to interfering amplitudes for the two paths that particles can take to reach 
the detectors. Intimately, bunching and anti-bunching reflect the statistics obeyed by the 
quantum particles: Bose- Einstein statistics imposes addition of the amplitudes, i.e., con- 
structive interference, while Fermi-Dirac statistics yields destructive interference. In the 
latter for elementary bosons, an anti-bunching effect signals that the matter state 

las no classical analog. Note that fermion anti-bunching was also observed for electrons 

9! [lol. H] and neutrons 12] . 

One encounters a more complex situation while studying composite bosons made of two 
fermions. In an ensemble of such quantum particles, the Pauli exclusion principle induces 
fermion exchanges between composite bosons which therefore do not exactly follow Bose- 
Einstein statistics. Semiconductor excitons constitute a good example of such composite 
bosons. These are made of electron-hole pairs and fermion exchanges between the electrons 

n 

and holes of these excitons are highlighted by optical nonlinearities [Ujj]. Interacting Fermi 
gases, made of e.g. 6 Li and 40 K, constitute other systems where composite bosons can be 



studied. In ultra-cold samples 



14J, these offer an interesting route towards precise inves- 



tigation of the Pauli exclusion principle. Indeed, exchanges between fermionic components 
can be varied, for instance when imbalanced spin populations interact [l5| |. 

Recently, a significant breakthrough was made in the theoretical description of composite 
bosons made of two fermions. A novel formalism, free from any mapping to an ideal boson 
subspace, has been constructed, and a visualization of the physical processes taking place 
between composite bosons (cobosons) has been proposed through the so-called "Shiva" dia- 



grams (see 



16| for a general review). This coboson many-body theory relies on two sets of 



2x2 scatterings: the "interaction scatterings" and the "Pauli scatterings" . The first ones 



3 



correspond to interactions between the fermionic components of two cobosons in the absence 
of fermion exchanges while Pauli scatterings correspond to fermion exchanges between two 
cobosons, without any fermion interaction. 

In this work, we use this new many-body approach to study the influence of fermion 
exchanges in statistical signatures of composite bosons. To highlight corresponding effects, 
we extend number-states to composite bosons by considering a many-body state constructed 
from N identical coboson creation operators. Hence, we evaluate the number of coincidences 

(2) 

in two-particle correlations, g N , as well as the variance of the field which is calculated 
through the Mandel parameter Qn- 

For number-states made of elementary bosons, \iPn) = Bq N \v), where Bq is the field- 
particle creation operator and \v) the vacuum state, it is known that gff and Qn read as (1- 
1/N) and (-1) respectively: Number-states exhibit an anti-bunched second-order correlation 
function and yield the greatest negative value for the Mandel parameter. This implies 
that measurements of such states are made with the greatest sensitivity, i.e., that noise is 
reduced at the smallest level. For composite bosons in a many-body state \ipN) — Bq\v), 

(2) 

we find that fermion exchanges modify the number of coincidences which now reads g N ~ 
1 + [—1 + 0(rf)]/N in the large N limit. 0(rj) is a positive correction whose dominant 
term increases linearly with the dimensionless parameter rj = N(cib/L) d associated to the 
composite boson density, a# being the coboson extension, L the sample size and D the space 
dimension. Consequently, two-particle correlations of composite bosons present a reduced 
anti-bunching effect when compared to elementary bosons. Furthermore, the amplitude of 
the Q-parameter for the state \i/jn) reads as Qn — —1 + 0(r] 2 ). Hence, the dispersion of the 
probability distribution is increased by fermion exchanges. In general, our analysis shows 
that it is necessary to take into account the Pauli exclusion principle between fermionic 
components of composite bosons in order to possibly extract statistical properties of these 
quantum particles in a precise way. 

I. PHYSICAL UNDERSTANDING 

Statistical properties of quantum fields are often studied through the second-order corre- 
lation function. The latter evaluates fluctuations in the number of field-particles for a given 
many-body state. Precisely, it measures the normalized probability to detect a particle 



i 



conditioned upon detection of a previous particle. The normalized number of two-particle 
coincidences in a state is defined as 

(2) _ {Bj 2 Bl) N 

where (A) N = (i^nIA^n) / (iPnIiPn) is the mean value of the operator A in state \i/jn), the 
field-particle creation operator being Bq. 

For elementary bosons, comparison between the normalized number of coincidences and 
1 allows us to deduce whether fluctuations of the number of particles follow classical or 
quantum statistics. Another way to study the variance of the number of particles and to 
characterize a quantum field is via the so-called Mandel parameter which precisely reads 

= (a% - m _ 

{n) N 

where n = BqB is the number-operator for the state 0. 

The parameter Qn compares the fluctuations of the field number-operator to that of a 
Poissonian source [4|. A negative Qn signals that the field statistics is sub-Poissonian and 
hence corresponds to the one of a quantum state. Most striking examples are obtained for 
the number-states |^jv) = Bq N \v) of particles with bosonic statistics, i.e. [B m ,Bj] = <S TOi . 
For these, we find gff = (1 — l/N) and Qjv — — 1 for any particle 0. Note that (—1) is 
the greatest possible negative value allowed for the Mandel parameter since the ratio in 
Eq.(2) is always positive: its denominator is the norm of Bq\i/jn) while its numerator is the 
norm of P±n\ij}N) where P± = 1 — \?Pn) (iPn\/ (iPn\iPn) is the projector over the subspace 
perpendicular to \^n), this operator being such that Pj_ = Pj_. 

To demonstrate that fermion exchanges play an important role in statistical signatures 
of composite bosons, we extend boson number-states to cobosons and consider the unnor- 
malized ket I^jv) = B^ N \v) to evaluate gffl and Qn. In order to understand how fermionic 
components affect these two quantities in a simple way, let us start with N=2 since most 
physical effects induced by the Pauli exclusion principle already appear with just two com- 
posite particles. 

A. Two composite bosons 

Since Bq\iP 2 ) is a zero-pair state, {^ 2 \Bq 2 Bq]^) is equal to (^ 2 |-6o \ v )( v \Bo\i^2)] so that 
by noting that (v\Bq\iP 2 ) is nothing but (^IV^), the normalized two-particle coincidences 
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reduce to 

(2)_ (jfe | fa) 3 / o\ 

92 ~ { ) 

To better see differences induced by the particle composite nature, let us brieffy reconsider 
elementary bosons, i.e., bosons such that [B m , Bj] = <5 mi . We then have B \ip2)=2\ijji) so 
that Sq|^2)= 2|u) and n\^p2) = 2 1 • Consequently ($2! $2) = 2 while (^l^fe) = 4 and 
(V^l^ 2 !^) = 8. This leads to the expected results, namely Q2 = —1 and g% =1/2 which 
is nothing but (1 — 1/N) taken for N = 2. 

We now turn to composite bosons. These have creation operators which do not exactly 



follow bosonic commutation rules 161 ] . but instead 

[B m , B\] — S mi - D mi , (4) 
where the so-called "deviation- from-boson" operator, D mi , is such that D mi \v) = while 

[D m ,B}] = J2(K n m ))+^j)Bl (5) 

n 

The parameters A(J^ are the composite bosons "Pauli scatterings"; they describe fermion 
exchanges between in the absence of fermion interaction (see Fig. l.(a)). 

The previous commutators give Bo\ip 2 ) = 2|^i)-2L2|f ) and B^L^ipi) = A2 1 ^1 ) + (L\ ~ 
2L\)\v) where we have set A2 = A([j [j) while 

4 = ^[A)o,4] = E A (oo°)4 (6) 

4 = |[A»,4]=£a(S£)A(?o°)4 (7) 

This allows one to find Bq\iI) 2 ) — 2(1 — A 2 )|t>) such that (^2 1 ^2) — 2(1 — A 2 ). This also leads 
to n\ip2) = 2^2) — 2L 2 |'0i). In order to calculate (^l^l^) i n an easy way, it is convenient 
to first note that (^l-^ll^i) = {' l l ) 2\B} i L\\v) . From (^ 2 |-Bo as gi ven above, we then find that 
(ip2 1 L\ I ^l) = 2(A 2 — A 3 ) where A 3 = Z) n A(° °)A(q °) corresponds to the Shiva diagram for 
fermion exchanges between three composite bosons 0, shown in Figure l.(b). All this leads 
to 

<^ 2 |^|^2> = 4(1 - 2A 2 + A 3 ). (8) 

The most convenient way to obtain (V^l^ 2 ^) is to first calculate BQh\ip2) knowing BqL\\iI)i) 
given above. By noting that (^l^l^i) — 2(A 3 — A 4 ), where A 4 corresponds to the Shiva 
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FIG. 1: Shiva diagrams representing fermion exchanges between composite bosons. Fermionic 
components are represented by solid and dashed lines, (a): Fermion exchange between two "in" 
cobosons in states leading to the formation of two "out" cobosons in states (n,m). This 

diagram represents the Pauli scattering A(^ *•). (b): Fermion exchanges between three composite 
bosons. Two possible ways to decompose fermion exchanges between three "in" cobosons and 
three "out" cobosons (to, 0,0) are represented, as appearing in Eqs (7) and (38). For m = 0, this 
corresponds to the three-body scattering A3 = J2 n ^(0 n)Mo 0) appearing in Eq.(8). 

diagram for fermion exchanges between four composite bosons 0, we end with 



Before going further, we can check that the above results for cobosons reduce to the ones 
for elementary bosons when all fermion exchanges, i.e. all A n 's, are dropped. 

For cobosons with center of mass momentum Q and relative motion index u , it is 
possible to show that A n = I (kj^o) | 2n - This makes A n a positive constant of the order 
of (a B /LY n -V D where as is the spatial extension of the relative motion state \vq) and D the 
space dimension. For 3D excitons or Hydrogen atoms, for which |(k|z/ )| 2 =647r(aB/-^) 3 /(l + 
k 2 a 2 B ) 2 , we find 



This leads to A 2 =(337r/2)(aB/L) 3 and A 3 =(41997r 2 /8)(aij/L) 6 which makes A 3 <A 2 and 
(^3-^2) appearing below, in the Mandel parameter, a positive constant equal to 
(202l7r 2 /8)(a B /L) 6 . 




(9) 




(10) 
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By collecting all the above results, we end with a normalized number for two-particle 
coincidences of composite bosons which reads as 

while the Mandel parameter is given by 

Q 2 ~-1 + (A 3 -A*). (12) 

This shows, for N = 2, that the number of two-particle coincidences is slightly higher 
for composite bosons than for elementary ones. We also find that the Mandel parameter is 
increased. As seen from Eqs.f lll|fT2l . the corresponding variations are solely controlled by 
the Pauli exclusion principle, through the exchange scatterings X n between n cobosons 0. 
Unlike elementary bosons for which the second-order correlation function and the Mandel 
parameter depend on N only, their counterparts for composite bosons also depend on the 
state at hand. 



B. Guess for N composite bosons 

If we now turn to N elementary bosons, number-states are such that u^n) = N\^n) 

(2) ~ 

such that g N = 1 — l/N and Qn = —1 for all N, whatever the state is. 

For composite bosons, a "rule of the thumb" gives intuitively the leading term in the 
small density limit, through the replacement of (as/L) D by (N — 1){(Ib/L) d in the result 
for N = 2. Since A3 and A| are both of the order of (as/L) D , we are led to guess the Mandel 
parameter as 

Qn — — 1 + (N — 1) 2 (A 3 - \\) ~ -1 + 0( V 2 ), (13) 

where rj is the previously defined dimensionless parameter associated to density. The situa- 
tion for gff is more ambiguous. From Eq. liTTj) . we might guess 

9$ » (l - ^) (1 + (N - 1)A 2 ) » (l - I) (1 + 0( V )) , (14) 

This would lead to a crossover from anti-bunching to bunching when l/N crosses 77, i.e., 
when the system size increases - which is very unlikely physically. However, we might as 
well guess 

9 ^ 1+ -H(;-D^ H ^ (15) 



s 



which maintains the anti-bunching effect for all sample sizes since the A n expansions per- 
formed in these calculations only hold for rj <C 1. The purpose of the next section is to 
demonstrate that Equations (TT3j) and (1T51) are indeed the correct expressions of Qn and g N . 
To do it, we must go deeper into the composite boson many body theory. 



II. SECOND-ORDER CORRELATIONS FOR N COMPOSITE BOSONS 



A. Key results from the coboson many-body theory 

Many-body effects between N composite bosons linked to fermion exchanges are obtained 



through the two following commutators 



[B m , Bj N ] = NB\ N ~\5 mi - D mi ) - N(N - l)^ 2 £ A(» , (16) 

n 

[D mu Bf] = NBf' 1 £ (AM + A(» J)) B\, (17) 



Note that these two commutators reduce to Eqs. (H|5l) when N = 1. 

Using them, it is possible to show that the normalization factor for N cobosons in the 
same state, differs from its elementary bosons value (N\), due to fermion exchanges taking 
place between them. This led us to write 

{v\B»Bl N \v) = MF N , (18) 

where the -Fat's, which enter all calculations involving a large number of identical cobosons, 
follow the recursion relation 



F N = Fjv-i — (N — 1)\ 2 F N _ 2 + (N - 1)(N - 2)A 3 F, 



jV-3 



vr nn-i (^-i)' A 



(19) 



The A n 's are the previously defined scatterings for fermion exchanges between N cobosons 
0, as shown in Fig. 1 for N=2 or 3. 
Equation ffT6|) readily shows that 

B \i> N ) = iV|^-i> - N(N - l)4|^iV-2>, (20) 
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with L\ given in Eq.(6). The above equation used for B \i[)n) and then for Ll\ipN-i), allows 
us to show that 

+ N - 1 

BlB^ N ) = \^ N ) + ^-^Bol^). (21) 

This last equation will turn very useful in calculating matrix elements involving the number- 
operator. 

B. Calculation of Qn 

Since (^n\Bq = (^tv+i|, Eqs. ffTH]) and (12TT) give the mean value of the particle number- 
operator as 

(h) N = i + ( N -l)^ = N+{N- 1)A« (22) 

where we have set = (F N+n — F N+n _i)/ F N . 

For N = 2 cobosons, Eq. (1191) gives F 2 = 1 — A2 and F 3 = 1 — 3A2 + 2A3, so that the above 
equation agrees with Eq. (jSJ). For N elementary bosons, F^ reduces to 1 and consequently 
all the A^ differences reduce to zero. Hence, the number operator mean value reduces to 
N, as expected. By contrast, Ajy for composite bosons are negative scalars since F^ is a 
decreasing function of N, as seen from Eq. (fT9]) . This makes the number operator mean value 
for cobosons smaller than its elementary boson value N, due to a "moth-eaten" effect similar 
to the one we have already found in other problems dealing with composite bosons: when a 
coboson is added to N other cobosons 0, the additional coboson feels the other N's through 
the Pauli exclusion principle. Therefore N elementary fermion pair states are blocked and 
thus missed in its linear combination, as if these were "eaten" by N "moths". This picture 
allows us to physically understand all decreases found for cobosons, when compared to their 
counterpart elementary boson values. 

By repeatedly using Eq. (l2"Tl) for (ip N \fi 2 \i/j N ) , split as ((^^\n)(n\ ; ip N )), we can show that 

(n 2 ) N = N 2 + (N 2 - l)Ag> + ^^A$. (23) 

Again, (n 2 ) ^ for cobosons is smaller than its N 2 value for elementary bosons, due to the 
same moth-eaten effect, since all the A^ are negative. Using F4 = 1 — 6A2 + 8A3 + 3A2 — 6A4, 
as deduced from Eq. ffl9l) . it is possible to check that the above equation agrees with Eq.Q 
when N = 2. 
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The two above equations allow us to write the mean quadratic deviation in a compact 
form in terms of two differences only, namely 

(n 2 ) N - {nf N — {N — 1) 2 [-Aff(l + A«) + j^^}- (24) 

Since, in the small density limit, A^ } « -iVA 2 + N(N - 1)(A 3 - Xf) and A^ } « -(TV + 
1)A 2 + {N + l)iV\3, as deduced from Eq. ffT9j) . we find that the mean quadratic deviation 
has no term in (a/L) D , its small density dominant term being N(N — 1) 2 (A 3 — A|). Using 
Eqs. (1221) . it is then easy to show that the small density value of Qm given in Eq. (1131) . as 
guessed from the rule of the thumb, is fully correct. 



(2) 



C. Calculation of g N 



The calculation of the normalized number of coincidences is a little more demanding. To 
evaluate the numerator of g N , we first rewrite BqBq using the commutator (j4]). This yields 



(Bl(BlB )Bl) N = (n 2 - n) N + (bId 00 B ) n . (25) 

The first term readily follows from Eqs. (122|23jl . To express the second term in a compact 
form is more difficult. Non-trivial manipulations are indeed necessary, otherwise we end 
with Fjv expansions which are far from obvious to sum. These manipulations are reported 
in the Appendix. They lead to 



(BlD 00 B ) 



N 




where R, defined as 



R 



1 A 



(2) 
N 



1 



iViV+1 A^(A^-1^ 
tends to \\ in the large sample limit (see Appendix). 
Using Eqs. (22,23,25,26), it is then easy to show that 

N 



(A 2 + A«) 



(26) 



(27) 



Bl 2 Bl) N = N(N - 1) 



A (i) 

N 



N+1 



-A 



(2) 
N 



2R 



(28) 



If we now remember that, in the large sample limit, A^ tends to —NX 2 , A^ to —(N + l)\ 2 
and R to A 2 ,, we find that the bracket in the above equation tends to [1 — (2iV — 3)A 2 ]. We 
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then deduce using Eq.(22) that g N 



is given by 



^(1-^)(1 + A 2 )^1 + 



-1 + NX 2 
N 



(29) 



in agreement with the result guessed in Eq. (1151) . since (N — 1) ~ N in the large N limit. 



By extending number-states to composite bosons, we have shown that the Pauli exclusion 
principle modifies statistical signatures of composite bosons many-body quantum states. 
This is underlined through the particular evaluation of the number of coincidences in two- 
particle correlations and the dispersion of the corresponding probability distribution. These 
explicitly reflect fermion exchanges between components of composite bosons. We find that 
the number of coincidences in second-order correlations is enhanced compared to elementary 
bosons, so is the variance of the field. However, main statistical signatures are found to be 
preserved. Indeed, the extension of number-states to composite bosons still presents a strong 
quantum character, with a Mandel Q-parameter close to (-1). 

To highlight these effects in the simplest way, we have first considered N = 2 composite 
bosons. From results obtained for just N = 2 cobosons, we can use a rather intuitive 
"rule of the thumb" to determine the number of coincidences in two-particle correlations as 
well as the Q-parameter for arbitrary N. The corresponding expressions, guessed from the 
N = 2 results, are thereafter confirmed by iV-body calculations. These are done following a 
procedure we recently proposed [3] to handle fermion exchanges between composite bosons, 
exactly. 

In general, our analysis confirms that the statistics of bosonic fields is modified by the un- 
derlying fermionic components of composite bosons. Precisely, fermion exchanges correlate 
composite boson states. This, in particular, makes the mean value of the number-operator 
smaller for composite bosons than for elementary ones. Therefore, coherence properties of 
composite bosons can not be directly deduced from their elementary boson counterparts. For 
that purpose, the derivation of operators peculiar to composite bosons, in order to possibly 
define "number-states" and "coherent-states", would be highly valuable. 



III. CONCLUSIONS 
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IV. APPENDIX 

This appendix is dedicated to the calculation of (BqDq B )n in a compact form. For 
that, we first use the commutator (jSJ). This allows us to write 

(BlD m B ) N = (D 00 BlB ) N - 2{L\B Q ) N . (30) 

The part (-^/vl-Doo i n the first term is obtained by mixing Eqs. (117j) and (120j) . This leads to 

D m \ip N ) = 2NLt\t/; N _ 1 ) = 2\t/; N ) - ^^—BoI^n+i). (31) 

By using Eq. ( 12TI) . it is then easy to show that the first term of Eq.( l30l) also reads as 

(DooBlB ) N = -2A« - 2 iV ^ V + ~ 1) Ag ) . (32) 

In the second term of Eq. fl30|) . we use Eq. (l20]) for BqI^n) and then the same Eq. fr20|) . but 
for L^I^at-i). This leads to 

(L\B ) N = -Aff - N(N - 1)R, (33) 

where R is defined by 

R = (iPn\l1 2 \iPn-2)n / 34 ^ 



All this allows us to rewrite Eq. (130l as 

(BlD 00 B ) N = - 2N jf + ~ - A^ + 2N(N - 1)R. (35) 

To get R in a compact form is considerably more difficult. Since it contains two Pauli 
scattering at least, through the two L\ operators, we can already say that its leading term 
must be in (a/L) 2D . This is going to give a rf contribution to gffl, negligible compared to the 

(2) 

dominant term of g N , expected to be in rj. Nevertheless, let us now show, for completeness, 
how a compact expression of R can be obtained. 

The trick is to note that, for L\ given in Eq.([n]), we have, using Eq. (1161) 

LMn) = £Mo r o) (nS^s-i) - N(N- l)EMn o°)^l^- 2 )) 

n \ m / 

= N\ 2 \i/>n-i) - N(N - 1)41^-2), (36) 

where L\ is nothing but the operator defined in Eq.([7]) since A(§ o)A(™ °) is the three-body 
exchange scattering with three cobosons on the right and two cobosons plus one coboson 
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m on the left, as readily seen from Figure l.(b). By using the above equation for (^1-^2 
but Eq.(j20l) for lI\iPn), it becomes easy to show that 

N-l F N + {M^n) ' 1 ' 

The next step is to calculate (^a^I-^sI^tv-i)- This is done by again using Eq. (l36|) for 
(^Ar_ 2 |L 3 and then Eq.( l20i) for L^jv-i)- This leads to 

(lp N -2\L S \^N-l) = j^j(^n-i\iPn-i) ~ jV2(jV _ i) (^IM + n 2 (N 2 - i^ ^ N+1 ^ N+1 ^ 

(38) 

So that we ultimately find the expression of R given in Eq. fl27|) . By using the values of 



and A^ in the small density limit, given above, we can show that R ~ X 2 ,. The correction 
to (BqD 00 Bq) n induced by this R term, is thus found to be of the order of t] 2 , as expected. 
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